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ABSTRACT A theoretical analysis of the ef- 
fect of force and torque on spontaneously twisted, 
fluctuating elastic ribbons is presented. We find 
that when a filament with a straight center line 
and a spontaneously twisted noncircular cross 
section is subjected to a sufficiently strong ex- 
tensional force, its average elongation exhibits an 
asymmetric response to large over and under- 
twist. We construct the stability diagram that 
describes the buckling transition of such ribbons 
under the opposing action of force and torque and 
show that all the predicted behaviors can be un- 
derstood in terms of continuous transformations 
between straight and spiral/helical states of the 
ribbon. The relation between our results and ex- 
perimental observations on DNA is discussed and 
a new reentrant spiral to rod transition is pre- 
dicted at intermediate values of twist rigidity and 
applied force. 



The first experiments on the stretching of double 
stranded DNA (dsDNA)§ were successfully described 
by the wormlike chain model |j| that accounts for the 
bending rigidity of the molecule . While an extension of 
this model reproduced the observed bell-shaped curves 
that characterize the response of dsDNA to torque in 
the weak stretching regime M, it could not capture the 
observed under/overtwist asymmetry at larger applied 
force [|[ E|. The feeling that linear elasticity alone can 
not account for the above observations, led to the devel- 
opment of a semi-microscopic theory in which the two 
strands of the double helix were modeled as elastic fil- 
aments connected by rigid rods that represent the base 
pairs, and an energy penalty was introduced for fold- 
ing the base pairs with respect to the central axis of the 
double helix[Q. This theory was able to reproduce most 
of the experimental observations on the deformation of 
dsDNA. However, the question is whether one is indeed 
forced to resort to such microscopic models to account for 
the experimental findings. In this paper we show that 
the linear theory of elasticity of spontaneously twisted 
ribbons combined with statistical mechanics, reproduces 
the qualitative features of the experiments on stretching 
and twisting of dsDNA, and provides important insights 
about the physical origin of the observed phenomena. 



1 THE MODEL 

The configuration of a thin ribbon of length I and asym- 
metric cross section is described by a triad of unit vec- 
tors {tj(s)}, where s (0 < s < I) is the contour dis- 
tance, t3 is the tangent to the center line and ti and 
t2 are oriented along the principal axes of the cross sec- 
tion. The orientation of the triad as one moves along 
the filament is given by the solution of the generalized 
Frenet equations that describe the rotation of the triad 
vectors, t-(s) = £\ fc e ljk LOj(s)t k (s), where e i]k is the 
antisymmetric tensor, the prime denotes differentiation 
with respect to s (t^ — dti/ds), and {u>i(s)} are the cur- 
vature and torsion parameters^. These equations can 
be rewritten in terms of the Euler angles 0, ip and ip: 

0' = lui sin ip + u>2 cos ip, 
ip' sin6> = — u>i cos ip + u>2 sin-0, (I) 
ip' sin 9 = (uji cos ip — u>2 sin ip) cos 9 + lu 3 sin 0. 

Notice that while the angles 9(s) and ip(s) describe the 
spatial conformation of the center line, ip(s) describes the 
rotation of the cross section about this center line. 

The elastic energy of a deformed ribbon is given by the 
sum of bending and twist contributions, U e i = Ubend + 
Utwist where [|| 
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Here k is the Boltzmann constant, T is the temperature, 
the bare persistence lengths ai and 02 represent the bend- 
ing rigidities with respect to the two principal axes of in- 
ertia of the cross section, and 03 is associated with twist 
rigidity. We assume that the stress-free reference state 
corresponds to a ribbon with a straight center line ori- 
ented along the x-axis and a cross section that is twisted 
about this line at a rate lu 3 q (this defines the twist number 
in Two = Iu)3o/2tt and the total angle of twist IuJzq). The 
above expression for the energy is based on linear theory 
of elasticity and applies to deformations whose charac- 
teristic length scale is much larger than the diameter of 
the filament |j) . Since we consider the deformation of the 
ribbon by forces applied to its ends, the total energy con- 
tains an additional term — (kTf/l) J ds sm9 cos tp where 



f 



/ is the force in units of kT/l. This theory is a gener- 
alization of the wormlike chain model to the case of a 
ribbon with noncircular cross section and nonvanishing 
spontaneous twist. 

In this work we only consider small deviations of the 
Euler angles 59, dip from their values 9 = 7r /2, fo = 
in the stress-free state. Physically, this corresponds to 
the strong force regime, / >■ 1, in which the tangent 
vector t 3 (s) undergoes only small fluctuations about the 
x axis. No restrictions on the magnitude of the devi- 
ation 5ip from the spontaneous value ^po(s) = ^30* are 
imposed. Expanding the energy to second order in the 
deviations 69, 6ip and introducing the complex variable 
$ = [66 + iSip] e-W yields 
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where we defined $ = <&' + iip'3>, and where UJ3 = 

<j)*<j> — <j><j>* 



We proceed to calculate the free energy of the de- 
formed ribbon, for a given value of the total angle of 
twist, ip(s) — ip'a, where we assume that rate of change 
tp' is constant (a priori, we do not exclude the possibility 
that solutions with different values of this constant coex- 
ist in the filament). The constrained partition function 
is evaluated by carrying out the functional integrals over 
$, $* with the Boltzmann weight exp(—U/kT). Since 
Utwist contains quartic terms in these fields, in order to 
calculate the integrals over <I> and <&* we introduce the 
Hubbard-Stratonovitch transformation of the twist con- 
tribution to the partition function 



exp (-Utwist/ kT) 



exp / ds 
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where 7 can be interpreted as a fluctuating torque. The 
resulting free energy is a quadratic form in <E> and and 
the Gaussian integrals over these fields can be carried 
out exactly. For convenience, we assume that the field 
<f>(s) obeys periodic boundary conditions, = $(0), 

and diagonalize the free energy by expanding $(s) in 
Fourier series, $(s) = ^2 n $ n e' 2 ™ s / L . A different choice 
of boundary conditions would affect our results only in 
the weak force region (in this work we will only consider 
the range / ^> 1). The integral over 7 is calculated by 
the steepest descent method. This yields the free energy 
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where we defined Lk 
f 2 + 2f[B(n 2 + Lk 2 )- 



- iP(1)/2tt and where h„(f,-f) = 
jLk] + (n 2 - Lk 2 )[A 2 (n 2 - Lk 2 ) 



+ 2B 7 Lfc-7 2 ] with A = {2tt) 2 ^/E 1 ~^/1, B = (27r) 2 (ai + 
a2)/2l and C — (2ir) 2 a 3 /l. A relation between Lk and 
the torque 7 is obtained by minimizing the free energy 
with respect to 7: 
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This relation has a simple geometrical meaning. The 
twist number Tw — (l/27r) J Q dsujs^s) can be related to 
the torque 7 by steepest descent evaluation of the inte- 
gral in Eq. @ that yields Tw = Tw + 7/C. Invert- 
ing Eq. (|l|) and assuming \S9\ <C 1, the twist number 

can be written as Tw — (l/27r) L ds [ip' — S9d (Sip) /ds 
For 1 89 1 <C 1 the writhe number can be expressed as 
Wr = (1/2tt) f Q dsS9d(8ip) /ds and thus Tw + Wr = 
Tp(l)/27t. Since the sum of twist and writhe numbers is 
the linking numbcr flo|| , we conclude that our definition 
of Lk coincides with the standard definition of the linking 
number. Inspection of Eq. (Q) shows that the first two 
terms on the rhs of this equation give the twist number 
and that the third term is the writhe number. 

A quantity that can be readily measured in experi- 
ments is the mean elongation of the filament (average 
end to end distance), (R) — —(l/kT)dJ-/df. The depen- 
dence of the elongation on the total angle of twist (2irLk) 
and on the extensional force is studied below. 



2 RESULTS 

Effect of torque on elongation. In Fig. [I] we plot 
the average elongation (R) vs. the linking number, 
for a spontaneously twisted ribbon (Lko = Twq = 10, 
A = C — 10, B = 50). At relatively low values of the 
force, we obtain a symmetric bell-shaped curve (curve 
a), in agreement with reference 0). Notice that fluctua- 
tions shift the peak of the curve to Lk max ~ 9.4 which is 
somewhat smaller than the spontaneous value. As / is 
increased, the dependence of the elongation on the angle 
of twist becomes progressively asymmetric (curve b) ; (R) 
decreases linearly with overtwist but is nearly indepen- 
dent of undertwist throughout the range < Lk < Lko. 
At yet higher values of /, the elongation becomes nearly 
independent of under or overtwist, in a broad range of 
linking numbers (curve c). 

In order to understand how the average conformation 
of the filament varies with the linking number, in the in- 
sert to Fig. ^ we plot the torque 7 and the twist number 
Tw vs. Lk, for a value of / that corresponds to curve 
b in this figure. Both 7 and Tw increase linearly with 
Lk in the range that corresponds to the flat portion of 
curve b, and approach constant values in the range where 
linear decrease of elongation with degree of over or un- 
dertwist is observed. In the range where Tw increases 
linearly with Lk (see insert), the writhe number remains 
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Figure 1: The average normalized elongation (R) jl is 
plotted vs the linking number Lk, for a ribbon with 
Lk Q = 10, A = 10, B = 50 and C = 10 (oi ~ 1, a 2 ~ 100 
and a3 = 10). The extensional force corresponding to 
the three curves are a) / = 50, b) / = 350 and c) / = 
1000. The dependence of the torque (twist) on the link- 
ing number for / = 350 is shown in the insert. 



small and the application of torque to the ends of a recti- 
linear filament results mainly in twist of its cross section 
about a straight center line. At large values of overtwist 
(ALk = Lk—Lk max > 0) for which Tw approaches a con- 
stant value, further increase of the linking number leads 
to the appearance of large mean writhe and the filament 
undergoes a transition to a spiral or a helical configura- 
tion. As will be shown in the following, this transition is 
related to buckling of elastic rods under torque^) (in the 
presence of thermal fluctuations the buckling instability 
is replaced by a continuous change of shape with increas- 
ing torque). Further increase of Lk increases both the 
amplitude (the radius) and the number of turns of the 
spiral and leads to progressive shortening of the average 
length of the filament, in agreement with the behavior 
observed in Fig. |l|b. Inspection of the insert shows that 
linear variation of Tw with Lk takes place over much 
broader range of undertwist than of overtwist, indicat- 
ing that throughout the plateau region in Fig. [j]b, the 
removal of spontaneous twist takes place mainly by un- 
twisting a straight filament, leaving (R) unaffected. The 
shallow minimum in the plateau region of curve b is re- 
lated to the existence of a reentrant spiral/helix to rod 
transition that will be discussed below. 

Force versus elongation. In Figs. ||a and ||b we 

plot / vs. (R) for different values of Lk correspond- 
ing to undertwist and overtwist, respectively. For small 
deviations from Lkg, the curves are symmetric under 
ALk <-> —ALk. For larger deviations a plateau-like 
region is observed at intermediate elongations, that is 
reached more rapidly for undertwist than for overtwist. 
Since for each set of parameters we find a unique solution 
ijj' of Eq. (j^), the above plateau is not associated with 
coexistence of two phases. Rather, the observation of a 



region of nearly constant log / over a large range of elon- 
gations means that in this region of parameters / varies 
slower than exponentially with (R) (e.g., as a power). 
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Figure 2: a) Plot of / vs. (R) /I, for increasing degrees 
of undertwist. The broken line corresponds to a sponta- 
neously twisted ribbon (Lk = Lkg), and curves 1 — 3 to 
ALk = —1.5, —3 and —10, respectively, b) Plot of / vs. 
(R) /I for increasing degrees of overtwist. The broken line 
corresponds to Lk = Lko, and curves 1 — 5 to ALk = 1.5, 
3, 7, 10 and 15, respectively. All elastic parameters are 
as in Fig. [l]. 

Buckling under torque. The simple physical pic- 
ture behind the above observations is intimately related 
to the well-known buckling transition under torque and, 
in order to gain intuition about the different parameter 
ranges, we examine the stability of the straight filament 
under the combined action of tension and torque. We 
substitute Eq. (Ji|) into exp(—U/kT) and look for the 
conditions under which the quadratic form in $ is no 
longer positive definite. In the absence of thermal fluc- 
tuations, a straight filament becomes unstable against 
buckling under torque at h n (f,j) = 0. The shape of the 
filament following buckling is given by a combination of 
Fourier modes $„ and where n is found by minimiz- 
ing h n (f, 7) with respect to n. In the presence of thermal 
fluctuations h n (f, 7) > everywhere and the buckling in- 
stability is replaced by a continuous deformation of the 
mean conformation of the filament from a straight line 
into a three-dimensional curve. In the limit of strong 
force and torque (the regime in which (R) decreases lin- 
early with Lk—see curve b in Fig. [j]), the resulting shape 
is dominated by a combination of modes with wave vec- 
tors 2Tr(Lk±n) Jl. The case n = corresponds to a 
simple helix with period l/Lk and the smaller axis of 
inertia of the cross section (direction of easy bending) 
rotating in the plane normal to its symmetry axis (see 
Fig. @). 

In Fig. |^ we plot the line of buckling transitions in the 
/— Lk plane for the same elastic parameters as in the 
preceding figures. The thick portions of the solid line cor- 
respond to instability at n = (rod to helix transition) 
and the thin portion of the line corresponds ton ^ (rod 
to spiral transition). The four inserts describe the aver- 
age configuration of the ribbon, at the locations shown 
by the corresponding arrows on the stability diagram. 
Insert 1 describes a ribbon with a cross section sponta- 
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Figure 3: The stability diagram in / vs Lk plane, for 
elastic ribbon with the same parameters as in Fig. [I] (the 
straight twisted ribbon configuration is stable against 
buckling above this line). The configurations of the fila- 
ment in different regions of the diagram indicated by the 
arrows are shown in inserts 1 — 4. The arrows denoted by 
a, b and c on the lefthand side of the figure correspond 
to the appropriate curves in Fig. El 

neously twisted about a straight center line (no applied 
torque), and corresponds to the stable region above the 
minimum at Lk — Lka. Insert 2 corresponds to the un- 
stable region to the right of this minimum (overtwist), 
in which buckling transforms the straight twisted ribbon 
into a helix, with the smaller moment of inertia of its 
cross section oriented normal to the axis of symmetry of 
the helix. The origin of the broad symmetric minimum at 
Lk = Lko is intuitively clear-small over/undertwist with 
respect to the equilibrium configuration destabilizes the 
filament against buckling and a larger extensional force 
is needed to maintain the straight state. The fact that in 
the absence of torque the instability appears to take place 
at / — > rather at some negative value of the stretching 
force corresponding to the Euler instability under load, is 
a consequence of our periodic boundary conditions; since 
the ends of the filament are not fixed, for / — > there is 
an instability against rigid rotation of the filament. The 
smaller minimum at Lk ~ appears only in the pres- 
ence of large spontaneous twist (Lk 3> 1) and bending 
asymmetry (ai, 0,3 <C CI2). Insert 3 shows a typical config- 
uration in the unstable region just below this minimum 
- buckling produces an untwisted ribbon bent along its 
easy axis. Since the stretching of this filament is opposed 
only by the smaller of the bending rigidities, it can be eas- 
ily stretched into a straight configuration. The presence 
of even small deviations from Lk = creates a nonplanar 
configuration shown in insert 4, the stretching of which 
invokes both the easy (ai) and the hard (02) bending 
axes, and requires a larger extensional force (hence the 
minimum) . 

The arrows a, b, c on the left hand side of Fig. || re- 
fer to the corresponding curves in Fig. [lj Recall that 
while thermal fluctuations were neglected in the deriva- 



tion of the stability diagram in Fig. 0, they were included 
in the calculations leading to Fig. ^l] (their main effect 
is to broaden the stability line into an extended tran- 
sition region between straight and spiral/helical states 
of the ribbon). Case (a) corresponds to the weak force 
regime in which strong fluctuations smear out the stabil- 
ity curve in the vicinity of the minimum at Lk = Lko, 
resulting in a symmetric bell-shaped curve (see Fig. |l|). 
Cases (b) and (c) correspond to the strong force regime. 
For small overtwists, all excess linking number goes into 
pure twist and the filament maintains (on the average) its 
straight configuration. At larger linking numbers beyond 
the fluctuation-broadened stability line, the filament de- 
velops positive writhe and undergoes a continuous rod- 
to-helix transition. Further increase of Lk increases the 
radius and number of turns of the helix, resulting in a 
linear decrease of (R). When the ribbon is subjected to 
undertwist, the filament remains straight and "unwinds" 
over a much larger range of undertwist compared to over- 
twist (because of the pronounced asymmetry of the sta- 
bility line), giving rise to the plateau regions in curves 
b and c, in Fig. [|. Note that there is an intermediate 
regime between cases (a) and (b), in which undertwist 
leads to reentrant behavior. For small degrees of under- 
twist, decreasing Lk initially unwinds the spontaneously 
twisted but straight filament and, as the stability line is 
crossed for the first time, the ribbon deforms into a he- 
lix with negative writhe and (R) decreases with increas- 
ing undertwist. As the stability line emanating from the 
Lk = minimum is approached, the helix undergoes a 
reentrant transition into an untwisted rod and (R) in- 
creases (a trace of this behavior is evident in curve b of 
Fig. |l|). Finally, at yet higher undertwists the stability 
line is crossed again and a transformation into a nega- 
tively twisted spiral takes place, accompanied by rapid 
decrease of (R) . 



3 DISCUSSION 

In this paper we studied the response of spontaneously 
twisted fluctuating elastic ribbons to externally applied 
torque and extensional force. The analysis is based on a 
combination of linear theory of elasticity and statistical 
mechanics, and the only input from an underlying mi- 
croscopic level of description is contained in the values 
of the elastic constants and of the rate of spontaneous 
twist. The agreement between our results and those of 
reference JtJ, is not surprising since the linear theory of 
elasticity is the long wavelength limit of any physically 
reasonable microscopic theory of solid behavior. Even 
though we made no effort to adjust our model parame- 
ters to fit dsDNA|Tl|], ° ur Figs, [l] and || contain most 
of the qualitative features of the experimental observa- 
tions on twisted and stretched DNA in the intermediate 
range of force, 1 — 70 pN|^, |(| . An exception to this state- 
ment is the predicted rapid decrease of elongation in the 
limit of large undertwist (Fig. nib) that was not observed 
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in experiment, possibly because the predicted decrease 
in elongation is preempted by a microscopic structural 
transition of dsDNA (alternatively, the decrease in elon- 
gation may occur at yet larger degrees of undertwist not 
reached in the experiments). 

All the features observed in Figs. [I] and || can be un- 
derstood in terms of a simple physical picture based on 
the stability diagram, Fig. ||, that describes the buck- 
ling instability under the opposing actions of torque and 
extcnsional force. When torque is applied to a straight 
spontaneously twisted ribbon, its response depends on 
the direction of the torque relative to that of spontaneous 
twist, and on the magnitude of / (assumed to be fixed 
during the process). For small degrees of overtwist, the 
filament remains straight and its cross section is twisted 
in excess of the spontaneous value. As the transition re- 
gion around the stability line is reached, the straight rib- 
bon deforms into a spiral curve or into a helix, depending 
on /. When the ribbon is subjected to undertwist it un- 
winds while maintaining a straight center line over much 
larger range than in the case of overtwist. Eventually, 
the transition region is reached and the straight ribbon 
deforms into a spiral or a helix. For the range of elastic 
parameters studied in this paper (ai <C a 3 <C 02) the sta- 
bility diagram contains an additional minimum at Lk ~ 
and, for small enough /, further undertwist leads to a spi- 
ral to rod transition (for forces below the critical value 
corresponding to this minimum the spiral transforms into 
a periodically bent planar configuration). At yet higher 
degrees of undertwist the straight filament deforms again 
into a helix. For larger values of / (above the value at 
the local maximum at Lk ~ 3 in Fig. ^J), the applica- 
tion of undertwist leads to the untwisting of a straight 
filament, resulting in a broad plateau in (R). Such be- 
havior was observed both in experiments on dsDNA[^[ |(| 
and in our Fig. [|. Note that although the value of / 
corresponding to curve b in Fig. ^ is well above the local 
maximum of the stability line of Fig. [|, the nonmono- 
tonic variation of (R) with degree of undertwist suggests 
that the fluctuations of the filament (and hence its elon- 
gation) are affected by the presence of the minimum at 
Lk ~ 0. While there is no compelling experimental evi- 
dence for reentrant behavior in dsDNA to date (however, 
a shallow minimum in (R) at intermediate values of /, 
is clearly visible in Fig. 3 of reference ||), it will be in- 
teresting to look for it in other systems such as dsDNA 
"dressed" by attached proteins, RNA, etc.. 

We would like to comment on the possible ramifica- 
tions of this work. Both the strength and the limitation 
of the present approach are its generality - while the 
description of the deformation of spontaneously twisted 
ribbons is a hitherto unsolved problem that is interest- 
ing in its own right, one may question whether it is a 
suitable model for dsDNA. Thus, one may ask whether 
the large asymmetry of the bending persistence lengths 
(0-2/0,1 ~ 100) necessary to produce the curves in Figs, 
^-^is physically reasonable. In order to answer this ques- 
tion notice that if we model a double helix by a sponta- 



neously twisted ribbon, the bending coefficients of this 
ribbon correspond to a hypothetical untwisted state of 
dsDNA in which two strands connected by base pairs 
form a ladder-like structure (two multiply connected, 
straight parallel lines). Such a structure is expected to 
have large asymmetry for bending in the plane of the 
ladder (02) or normal to it (ai) (dsDNA was modelled as 
a ribbon polymer made of two semi-flexible chains that 
could not be bent in the plane of the ribbon, in reference 
PH). The effective persistence length a e ff of our model 
is given by the relation a~h = (a^ 1 + a 2 1 )/2 and ap- 
proaches 2cii in the limit 02 S> a-\ (thus, while a± ~ 25 nm 
is determined by the persistence length of dsDNA, 02 has 
no direct physical interpretation). 

Another issue involves the description of deformation- 
induced internal transitions in dsDNA invoked by 
experimenters |L3| in order to explain the observed 
plateau in plots of / vs. (R), as well as the observed 
opening of the DNA bases at large undertwists (by hy- 
bridization with complementary fragments), and the in- 
creased reactivity of some bases at large overtwists. Even 
though in our continuum elastic model, there is no direct 
reference to internal structure, analogies between shape 
changes and structural transitions can be readily drawn. 
Thus, variations of internal structure would correspond 
to changes of the average conformation of the filament 
and it is plausible that the conformations shown (see in- 
serts 1 — 4) in Fig. || differ in their chemical reactivity, 
consistent with the observations of reference Fur- 
thermore, the experimental observation that even rela- 
tively small undertwists lead to denaturationjll|, is con- 
sistent with our result that over a large range of under- 
twists, the cross section unwinds around a straight center 
line and no writhe develops (one expects twist around a 
center line to be more effective than writhe of the center 
line in disrupting the interaction between the base pairs) . 
A deeper difference between our discussion of force and 
torque induced shape changes and that of transitions be- 
tween different internal states of dsDNA in reference E| 
is that while we find only continuous transformations of 
average conformation, discontinuous first order internal 
transitions (from B to P and S forms) and phase coexis- 
tence are invoked in the latter work. 

Since our model assumes small deviations of the Eu- 
ler angles from their stress-free values, it can not de- 
scribe large amplitude defects (with SO, Sip > w) such 
as plectonemes. However, our analysis shows that for 
/ 3> 1, the formation of plectonemes is preempted by 
the appearance of homogeneous spiral and helical struc- 
tures (recall that the classical torsional buckling insta- 
bility corresponds to a transition to a helical state with 
arbitrarily small amplitude [Qj), even though plectonemes 
may appear at yet larger applied torques. Finally, the 
fundamental limitation of our theory is that since it is 
based on the linear theory of elasticity of thin inextensi- 
ble rods, it can not describe the observed overstretching 

and one 
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of B-DNA by 1.7 times its native length) 14, 
must resort to microscopic fq or thermodynamic Jl6| ap- 
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proaches. 
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